V 


I 




* 


V 


tv. 

* 

.-jT 


4$ 


Mr 0, 


. ^ 


£ 



y, 





CEPilicijjj 

Pclytec hni c“Lto 0fl e \ IO p G “" U , ! 7 _'^-«l«r 

CSCL 07D 




V 



63/06 


.«J 


N74-22805 


Unclas 

39401 

I 


R.P.I. Technical Report MP-4l 


COMPOSITION DEPENDENT EFFECTS IN 
GAS CHROMATOGRAPHY 

Raymond C. Lavoie 

National Aeronautics and Space 
Admini stration 


Grant NGL 33-018-091 


Analysis and Design of a Capsule Landing System 
and Surface Vehicle Control System for 
Mars Exploration 


Rensselaer Polytechnic Institute 
Troy, New York 


May 1974 



1 


ABSTRACT 

Continuing research sponsored by NASA is concerned with 
developing fundamental concepts which are required to optimize 
a gas chromatograph - mas6 spectrometer chemical analysis system 
suitable for use on an unmanned roving vehicle for Mars explor- 
ation. Prior efforts have developed simulation models for the 
chromatograph which have been compared with data obtained from a 
test facility. Representation of binary systems by superposition 
was previously shown to be a first-order approximation and in 
certain cases large discrepancies were noted. Thl6 subtask 
had as its objective the generation of additional binary data and 
analysis of the observed nonlinear effects. 

Modifications were completed on the test facility and data 
analysis programs to improve accuracy of the results. Data were 
gathered on the porous C-102 column and the nonporous C-1500 
column for pentane, heptane, acetone, methyl ethyl ketone, and 
ethylene at various temperatures and system conditions. 

The thermodynamic parameter mRo for heptane on C-102 was 
found to vary linearly with component concentration in the gas. 
This paralleled the nonlinearity observed in the pentane-heptane 
binary system on C-102. Simulation was improved as the system 
composition became more dilute. A similar effect was observed for 
ethylene on C-102. Experiments with the pentane-heptane on a 
DES column showed the system to be well modeled by the equilibrium 
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adsorption model and superposition. This demonstrated that the 
extent of the composition dependent effect was specific to each 
chemical-adsorbent system. 

The theoretical Langmuir isotherm was examined as a possible 
means to account for this composition dependent effect. Orthogonal 
collocation was used to solve the equilibrium adsorption model 
using the Langmuir Isotherm instead of the linear Isotherm. 

The resulting solutions show qualitatively that the isotherm 
models actual column behavior. Quantitatively, simulation was 
inadequate and further experimentation with the numerical solution 
of this nonlinear partial differential equation is required. 
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I. INTRODUCTION 

' 

This report describes recent work done on multicomponent 
gas chromatography as part o£ an effort to design chemical 
analysis equipment for the unmanned martian roving vehicle 
being developed at Rensselaer* The overall task is to provide 
design criteria and engineering techniques to build an optimum 
system. The plan 16 to couple a gas chromatograph used to 
separate mixtures, with a mass spectrometer capable of analyzing 
chemical components. This system must be flexible enough to 
perform a variety of tasks including atmosphere analysis, pyro- 
lytic tests, and incubation tests of soil and atmosphere samples 
to detect possible biological activity on mars. The most crucial 
design specification for this system is that it be small enough 
and light enough to fit on board the vehicle. 

Gas chromatography is a means of separating a mixture of 
differing ch mical species by utilizing the dynamic mechanism . 
of adsorption - desorption between a packed bed and a flowing 
inert gas stream. Separation occurs because of differences in 
the adsorption isotherms of the various components. The more 
strongly a component is adsorbed the longer its elution time 
will be. In an ideal system an input pulse containing a mixture 
of components will elute as a series of pulses each containing 
one of the components. 
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A gas chromatograph can be used for much more than separ- 
ating mixtures however. A composition detector at the outlet 
of the column can monitor the area of each component pulse. 

This area is related to thi percentage of that component in the 
original sample. The relative amounts of the components can 
be calculated and a good guess at the identification of each 
species can be made, since each substance elutes from a given 
column at a characteristic time. The mass spectrometer provides 
positive identification. 

Since gas chromatography i« a complex process with many 
dynamic mechanisms continually taking place, mathematical models 
are being developed to predict system behavior under a wide 
variety of conditions. These predictions are compared to actual 
data. Deviations from experimental results point out the flaws 
in the models. Theoretically an exact model which would in- 
corporate all phenomena observed in chromatography could be 
developed. However, the complexity of such a model would make 
its numerical solution impractical even with present day computers. 

One must be satisfied with models describing only the 
most significant phenomena affecting the behavior of the column. 
These simple models can approximate reality to the extent that 
they are based on valid assumptions. If an important mechanism 
is not Included in the model, predicted results will deviate 
from experimental results in a consistent and recognizable 
way. The model must then be amended. In this way, the actual 
conditions are identified where each physical phenomenon is 



significant to the overall process* Simplified but valid models 
are the result* 

Special problems arise in the study of multicomponent chrom- 
atography. ?he possibility of interaction between the various 
components* and other nonlinearities must be considered. 1'he 
first approach to this problem is to assume no interaction 
between the components end to use linear superposition of in- 
dividually modeled responses to predict the overall response. 

In fact when this approach was tried large deviations from 
experimental data were found for one column* the chromosorb 102* 
The work presented here clarifies the cause of the deviations 
and explores the use of a model incorporating the nonlinear 
Langmuir isotherm to explain them. 



II. SUMMARY 


The gas chromatograph facility has been modified to improve 
the usefulness and accuracy of the data. A sleeve added to 
the liquid injection block has caused syringe injected input 
pulses to the column to be narrower than had been attained 
previously. A change in the circuitry of the detector controllers 
has made possible continuous attenuation of the output signal 
to the recorder. 

The computer programs used to analyze data have been updated 
to provide more information. The General Data Reduction Frogr^ 
was modified to report compositions as actual mole fractions 
and a moment analysis has been performed on the input and ou. t 
pulses for comparison wi th moments predicted from models. A 
program for modeling binary systems, using linear superposition 
has been written and can be used to predict the separation of 
two components for any set of conditions. 

New data have been gathered on the porous Chromosorb 102 
column and the nonporous Carbowax 1500. Binary samples could 
not be separated on the 0*1500 column. Experiments on a new 
nonporous column with a coating of di-2-ethylhexyl sebacate (DES) 
have shown that its behavior is extremely well modeled by the 
equilibrium adsorption model. Pentane-heptane mixtures can be 
separated on this column and linear superposition of pure com- 
ponent data models this behavior quite well. ► * 
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Experiments using various volumes of heptane on the Chromo- 
sorb 102 column have shown that the parameter mRo for this 
system is a function of the amount of component present in the 
inert gas stream and is not rongly affected by the presence 
of other components. The nonlinear Langmuir isotherm has been 
proposed for the model to account for this effect. 

This new model is nonlinear and cannot be solved analytically. 
The method of orthogonal collocation as developed by Woodrow 
was adapted to this problem. Results indie ‘e that the method 
in its present form is not suitable for quantitative analysis of 
the model. On a qualitative basis the solutions obtained indicate 
that a nonlinear Isotherm of the Langmuir type can predict the 
deviations from simpler models, observed on the C-102 column. 
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III. BACKGROUND 

A substantial amount of experimental work has been com- 
pleted on the chromatograph teat facility, over a wide range 
of temperatures, pressures, flow rates and input pulses, using 
several columns and chemicals with different physical rnd chemical 
characteristics. In this way a variety of situations has been 
studied where different physical phenomena gain or lose influ- 
ence on chromatograph behavior. So many phenomena are possible 
that only a complete model would be able to predict the observed 
results in every situation. 

• \ , previous investigators have contributed to the develop- 
ment of mathematical models to predict the behavior o.. gas 
chromatograph systems (Keba and Woodrow, 1972; Keisch, 1973; 

Sliva, 1968; Taylor, 1970; Voytus, 1969). The most useful of 
these models has been the equilibrium adsorption model which 
successfully predicts first order and in some cases second order 
behavior. This model takes into consideration axial interparticle 
diffusion, convection and equilibrium adsorption-desorption 
using a linear isotherm. The derivation of this model has been 
shown by Voy tus(1969). In its mathematical form, the model 
is represented by the following partial differential equation 
which is a mass balance for the Injected species in the fluid 
and solid phases of the column. 


(l + 1/mRo) <)y/<)o - 


d y/d z 
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Boundary Conditions: 

y(O,0) b input pulse 
lira y(z,0) =» finite 

The final solution of this equation is the following convolution 
integral: ^ 

y (2 ,0) » Jy(0, O-'H . y^.('T') d'f- 
vhere: ° 

y(0,G) c the input pulse 

y^ O") « the unit impulse response at the end of the 
column 

; - 1/2 JlW ex P 4-^“ } 

where: 

jQ - +■ l/i n<R-o) 

The Peclet number is a measure of the dispersion of the sample 
component due to axial diffusion in the column. Smaller values 
of Pe indicate greater dispersion. The thermodynamic parameter 
mRo determines the elution time of the Injected component. 

Larger values of mRo indicate that a component is not strongly 
adsorbed in the column and that it will elute at a sooner time. 
This model does not consider intraparticle phenomena, adsorption- 
desorption kinetics or nonlinear isotherms. The Importance of 
intraparticle phenomena is not yet well '..town. A study of the 
nonequilibrium adsorption model has shown that mass transport 
between the carrier gas and the adsorbent is too fast to have 
any limiting effects on systems studied thus far. The question 
of nonlinear Isotherms is discussed in this paper. 



The equilibrium adsorption model was solved using numerical 
convolution with the computations being performed with the General 
Data Reduction Program (GDRP). In the output of this program, 
model predictions are compared directly to actual data. In 
some situations, the model simulations predicted system behavior 
quite well, Indicating that all major dynamic phenomena had 
been accounted for. In a few situations the model predictions 
deviated consistently from actual data. 

The most notable deviation occurred with ethylene on chromo- 
sorb 102. Figure 1 compares the model to actual data and there 
1 8 little resemblance between the two. The actual output is 
very spread with much tailing, while the modeled response is a 
relatively sharp peak. When heptane was v ed on the C-102 
column there was also more spreading in the actual data than in 
the simulation. ! 

Another type of deviation was found when a binary mixture of 
pentane and heptane was used as the sample input on th& C-1Q2 
column. The tro species were separated but Meisch (1973) found 
that the parameter mRo was a significant function of composition. 
Figure 2 snows how the composition dependent mRo shifted the 
peaks of the actual data away from the mode? which is based on 
the linear superposition of pure component simulations. The 
figure also shows extra spreading of actual data that is not 
predicted by the model. 

The equilibrium adsorption model obviously does not represent 
all o i the physical behavior that is occurring in the column. 















There ere many possibilities to consider* Melsch's findings 
concerning the porosity of the column packing, support the notion 
that lntrapartlcle diffusion Is present. Woodrow (1974) is 
working on this question* The variation of mRo with composition 
supports the idea that a nonlinear Isotherm should be used in 
the model* This Idea Is Investigated In this paper* Perhaps 
both situations are Involved* The question is difficult to 
resolve because of problems encountered in solving the resulting 


mathematical models 



IV. RESULTS AND DISCUSSION 


A. Modification of Equipment 

All experimental gas chromatograph work has been done on 
a test facility which was described by Baer and Benoit (1971). 

Though this equipment had been proven adequate it needed improve- 
ment to increase the accuracy of the data it provided. Meisch (1973) 
pointed out that the ten position attenuator on the detector 
control was not a versatile enough adjustment to give full scale 
readings on the oscillograph recorder for the range of samples 
studied.. This limitation decreased the relative accuracy of the 
data. In most cases a setting of one caused the readings to 
go off scale while the next setting of two caused deflections 
less than half the width of the strip chart. The 'ten position 
Attenuator was replaced with a potentiometer, making continuous 
attenuation possible. The new circuitry is shown in Figure 3. 

In separating a mixture of components on a gas chromato- 
graph, it is desirable to use as narrow an input pulse as possible. 
This fact can be illustrated using moment analysis. It has 
been shown that the variance, or second moment about the mean, 
of the eluted peak is equal to the sum of the variances of the 
input pulse and the impulse response predicted by the model 
(Voytus, 1969). The width of a pulse is related to the square 
root of the variance. In order to ensure that a sample vlll be 
separated, the variance of the eluting peaks should be email 
to prevent overlapping. 




CARLE DETECTOR CONTROL CIRCUIT — FIG. 







A study of the liquid injection system showed; that the 

» 

peak of the input pulse occurred as much as 15 seconds after 

1 

the actual syringe injection. The cause of this spreading was 
found to be an excessively large vaporisation chamber in the 
injection block. As shown in Figure 4, the Injection chamber 
is formed by a cylindrical hole drilled from the top of the 

block to within a half inch of the bottom. A cylindrical heater 
is placed in the chamber, sealing the top of the block. It is 
immersed in the stream of helium which enters at the top of the 
block and leaves at the bottom. The sample is Injected through 
a septum which is located near the entrance of the inert stream. 
In order to reduce the dead volume a tubular sleeve was designed 

which fills most of the space between the heater and the chamber 

♦ 

wall. As is shown in Figure 4, the sleeve has slots which 
channel the gas stream directly through the chamber. To be 
sure that the liquid sample would completely vaporize in the 
smaller residence time, the wiring was changed to give continuous 
full power to the block heater. 

Figure 5 shows how the addition of the tubular sleeve to 
the Injection chamber has narrowed the input pulse. Under the 
same flow conditions with two equal sample sizes, the peak of 
the input peak occurs at 10 seconds before the change and at 
3 seconds after the change. There is also a substantial decrease 
in the variance which was 131. before and 25. after the sleeve 
was added. It should be noted that the characteristics of the 
input depend on many factors including flow rate, sample size 
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amount of time taken to depreaa the syringe, positioning of the 
syringe and condition of the rubber septum. It has been found 
for example that narrower input pulses result when new septa 
are used. 

B. Linearity of the thermistor probe 

While the accuracy of the data gathered on the gas chromato- 
graph facility depends on the equipment, the validity of these 
data depends on correct Interpretation of the information received. 
Data processing methods assumed linearity between the electrical 
signal produced by the thermistor probes of the detectors and the 
composition of the sample in the carrier gas. A simple experiment 
was devised to verify this assumption. 

The column was operated at constant flow rate and temperature. 
The electrical controls were set so that a 3.0 /* 1. sample of 
heptane would register a full scale deflection on 'the oscillo- 
graph recorder. This input pulse was recorded and subsequently 
a 1.5 jx 1, sample was injected and its pulse was recorded. A 
similar procedure was followed for a 2.0 /*1. and 1.0 /jl. sample. 
Linearity of the detection system would cause the area under the 
curves to be proportional to the sample size. The areas .were 
computed and the results follow in Table I. 

TABLE Z 

Data Concerning the Output of the Thermistor Probe 

Sample 1 Sample 2 Ratio of Areas 

3.0/jl 

2.0/a/ 


Test 1 
Test 2 


1.5/*/ 

1.0/W 


1.948 

2.073 
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The ratio of area* ahould be 2.0 In both caeea if the sample 
sizes are accurate and the detector Is linear* The actual results 
are about 3*0% low In the first run and 3.5X high In the second* 
Considering the limited accuracy in experimental technique and 
In computations of the areas, the deviations are not significant. 
The thermistor detector appears to be linear in the composition 
range of interest. 

C. Modification of Data Processing 

Some of the computer programs used for data analysis have 
been changed and new programs have been written to increase the 
amount of Information they provide* The Ceneral Data Reduction 
Program has been modified to read in the sample size of the 
input pulse and the flow rate of the carrier gas so that it 
now reports input and output compositions as actual mole fractions 
rather than as normalized concentrations* This information is 
useful in comparing compositions of input pulses of different 
sizes and shapes* It also can provide information concerning 
the concentrations of sample in the carrier gas at which com* 
position dependent nonllnearltles might be significant* Sub- 
routine STATMO, written by P* Woodrow (Appendix A), was added 
to the program. It calculates the first moment about the origin 
and the second and third moments about the mean for the input and 
output pulses, and reports them in both dimensionless and real 
time domains* These moments are useful when the moments of a 
model can be calculated from the Laplace transform of the model* 

By comparing the moments of the model to the actual moments , 
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a model can be quickly verified without obtaining a tirr,~ uomaln 
solution. This technique has proved feasible for a complicated 
lnter*intrapartlcle model (Woodrow, 1973). 

The method of analyzing binary data has been revised some* 
what* Figure 6 Is a schematic of the new procedure* The Binary 
Analysis Program, written by P. Keba (1972), requires the input 
of simulated and actual data points* "his program Itself has 
been changed only so that It reports compositions as actual 
mole fractions* A new program, SIM has been written which uses 
the Input pulse and system parameters of the binary mixture 
and the mRc of the pure components at the same temperature to 
simulate the system and punch the simulated chromatogram data 
points on data cards* The linear superposition of the pure 
component simulations Is compared with the actual data by the 
Binary Analysis Program* * 

A program has been written for simulating the separation 
of a binary sample using Information from pure component data* 
Figure 7 shows the organization of the program* System parameters 
such as column and packing dimensions as well as flow rates and 
types of Input pulses can be varied to find the effect of each 
parameter on the output* This program has found several uses* 

It can predict if a column can be used to separate a given 
mixture with the limitations Imposed by the experimental equipment* 
This information has guided the recent purchase of a new column. 

A revised version of this program has been used by Stovall (1974) 
to study resolution of partially separated peaks* 
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Figure 6 Procedure For Analysis of Binary Data 
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D. Summary of Data 

For convenience, all of the chromatographic data obtained 
to the present has been summarized In tabular form. The system 
parameters, temperature, flow rate, gas velocity, Peclet number, 

and mRo are reported for each run. In addition, the sample 
size and the time at which the input pulse reached a peak have 
been reported where possible to give a measure of the relative 
characteristics of different inputs. 

Table II shows data for pentane, heptane, ethylene, methyl 
ethyl ketone and acetone on chromosorb 102. New data on this 
column include heptane at varioue sample sizes and methylethyl 
ketone. Table III reports binary data obtained on C-102. 

Figure 8 illustrates the variation of m^o with temperature for 
mek on this column. 

Data for the carbowax 1500 column are summarized in Table IV. 
New data include acetone and mek run over a range of temperatures, 
and ethylene at many temperatures, flow rates, sample sizes and 
with syringe Injection. The variation of mRo with temperature 
for pentane, heptane ethylene, mek and acetone are shown in 
Figures 9 through 13. 

The new DES column has been tested with pentane and heptane. 
Pure component data are given in Table V, and binary data are 
in Table VI. The effects of temperature on mRo for pe.itane 
and heptane are shown in Figures 14 and 15. The effects of 
sample composition in a binary mixture, and sample size on mRo 
are shown -n Figures 16 and 17. Dimensionless r ments for 
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TABLE II 

PURE COMPONENT DATA 
ON C11ROMOSORB 102 


PENTANE (150°, 175°, 200°) 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

N 

tog 

mRo 

°c 


sec 

cc/min 

cm/ sec 




150° 

2.0 

12.5 

29.4 

21.13 

9128 

82640 

0.0753 

175' 

2.0 

12.0 

30.0 

22.84 

8978 

83960 

0.1017 

200° 

2.0 

12.0 

26.4 

21.22 

7735 

96560 

0.1494 


HEPTANE <150°, 175°, 200°) 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

■ N. 

tog 

mRo 

°C 


sec 

cc/min 

cm/sec 




150° 

2.0 

11.0 

29.4 

21.13 

10690 

70200 

0.0177 

1-0° 

2.0 

9.0 

41.6 

25.00 

13770 

53690 

0.0184 

175° 

2.0 

8.0 

36.6 

23.30 

12140 

61640 

0.0383 

175° 

2.0 

12.0 

30.0 

22.84 

10520 

71340 

0.0320 

200° 

2.0 

15.0 

20.4 

16.40 

1 

102100. 

0.0711 

200° 

2.0 

11.0 

34.0 

22.85 


67570 

0.0698 









TABLE II (cont.) 


j 

1 

HEPTANE 200° C Varying Sample Concentration 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

N. 

tog 

mRo 

°C 

V 1 

sec 

cc/min 

cm/ sec 




200° 

0.2 

3.0 

28.6 

22.99 

9781 

76600 

0.0482 

200° 

0.5 

4.0 

28.6 

22.99 

9781 

76600 

0.0514 

200° 

1.0 

3.0 

28.6 

22.99 

9781 

76600 

0.0547 

200° 

1.5 

3.0 

28.6 

22.99 

9781 

76600 

0.0590 

200° 

2.0 

4.0 

28.6 

22.99 

9781 

76600 

0.0611 

200° 

3.0 

4.0 

28.6 

22.99 

9781 

76600 

0.0666 


ETHYLENE (20°, 28°, 50°, 75°, 100°, 125°, 150 , 175 ) 


Temp 

°C 1 

Sample 

Size 

Flow 

Rate 

cc/mln 

Gas 

Velocity 
cm/ sec 

Peclet 

No. 

N. 

tog 

mRo 

o 

o 

<M 

100 

43.0 

23.56 

9250 

83780 

0.1225 

28° 

100 

43.0 

25.65 

10210 

76270 

0.1117 

50° 

100 

43.0 

26.46 

9744 

79750 

0.1942 

75° 

100 

43.0 

27.79 

9294 

83430 

0.3199 

100° 

100 

43.0 

27.26 

8821 

87650 

0*4230 

125° 

100 

43.0 

27.74 

8446 

91270 

0.4785 

150° 

100 

43.0 

28.33 

8102 

94870 

0.6008 

175° 

100 

43.0 

28.76 

7793 

98340 

0.5869 
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TABLE II (cont) 
Methyl Ethyl Ketone (150°, 175°, 200°) 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

N . 

tog 

mRo 

°c 


sec 

cc/min 

cm/ sec 




150° 

0.5 

1.0 

40.4 

26.46 

10850 

70240 

0.0371 

175° 

0.5 

.75 

40.4 

28.02 

10500 

72610 

0.0565 

200° 

0.5 

2.0 

42.6 

29.60 

10580 

72070 

0.0906 


Acetone (100° 

» 125°, 

150°, 175°, 200°) 




Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

N 

tog 

mRo 

°C 


sec 

cc/min 

cm/sec 




100° 

2.0 

12.0 

43.0 

27.26 

8689 

88960 

0.0286 

125° 

2.0 

1.75 

24.0 

19.84 

4938 

148700 

0,0576 

125° 

2.0 

12.0 

43.0 

27.74 

8309 

92760 

0.0527 

150° 

2.0 

12.0 

43.0 

28.33 

7965 

96480 

0.0968 

175° 

2.0 

12.0 

43.0 

28.76 

7652 

100100 

0.1219 

200° 

2.0 

12.0 

43.0 

29.32 

7388 

103400 

0.1734 


TABLE III 


BINARY DATA ON CHROMOSORB 102 


Pentane - Heptane 150°C 



Flow Rate 
cc/min 

Gaa Velocity 
cm/sec 

Peclet No. 

N 

tog 

Pentane 

29.4 

21.13 

9128 

82640 

Heptane 

29.4 

21.13 

10690 

70200 


Wt Fraction 
Pentane 

Wt Fraction 
Heptane 

Input 

Peak 

Pentane 

mRo 

Heptane 

mRo 

.99 

.01 

13.0 

0.0747 

0.0116 

.75 

.25 

13.0 

0.0725 

0.0137 

.50 

.50 

10.5 

0.0673 

0.0153 

.25 

.75 

11.0 

0.0663 

0.0163 

.01 

.99 

11.0 

0.0599 

0.0171 
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TABLE HI (cont) 


Pentane - Heptane 175° C 



Flow Rate 
cc/min 

Gas Velocity 
cm/ sec 

Peclet No. 

N„ 

tog 

Pentane 

30.0 

22.84 

8978 

83960 

Heptane 

30.0 

22.84 

10520 

71340 


Wt Fraction 
Pentane 

Wt Fraction 
Heptane 

Input 

Peak 

Pentane 

inRo 

Heptane 

mRo 

.99 

.01 

11. 0 

0.1022 

0.0243 

.75 

.25 

12.0 

0.0988 

0.0261 S 

.50 

.50 

15.0 

0.0945 

0.0275 j 

.25 

.75 

11.5 

0.0925 

0.0277 j 

.01 

.99 

11.5 

0.0895 

0.0295 i 


Pentane - Heptane 200° C 


- . -- 

Flow Rate 
cc/min 

Gas Velocity 
cm/sec t 

Peclet No. 

N tog | 

i 

i 

Pentane 

26.4 

21 

.22 

7735 

96560 ‘ 

Heptane 

26.4 

21 

.22 

9118 

81960 

Wt Fraction 

Wt Fraction 

Input 

Pentane 

Heptane 

Pentane 

Heptane 


Peak 

mKo 

mRo 

.50 

.50 


11.0 

0.1471 

0.0531 , 

.25 

.75 


12.0 

0.1472 

0.0553 i 

j 

.01 

.99 


11.5 

0.1398 

0.0582 | 


* 


I 







TABLE IV 


PURE COMPONENT DATA 
ON CAR BOW AX 1500 


ACETONE (75°, 100°, 125°, 150°) 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

N 

tog 

mRo 

°C 


sec 

cc/min 

cm/sec 




75° 

0.5 

3.0 

34.2 

27.57 

7654 

100000 

0.1315 

100° 

2.0 

10.0 

31.8 

24.90 

6757 

112100 

0.1875 

125° 

2.0 

10.5 

29.6 

24.80 

6051 

123900 

0.2342 

150° 

2.0 

10.5 

27.7 

24.6 

5457 

136000 

0.2632 


PENTANE (75°, 100°, 125°, 150°) 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

Mo. 

tog 

mRo 

°C 


see 

cc/min 

cm/ sec 




75° 

0.5 

3.5 

34.2 

27.57 

12020 

62930 

0.3333 

100° 

2.0 

10.5 

31.8 

24.94 

10780 

70280 

0.4566 

125° 

2.0 


29.6 

24.77 

9763 

77470 

0.4280 

150° 

2.0 

9.5 

27.7 

24.63 

8878 

84840' 

0.4541 
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TABLE IV (cont) 

HEPTANE (75°, 100°, 125°, 150°) 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gae 

Velocity 

Peclet 

No. 

tog 

mRo 

°C 


sec 

cc/mln 

cm/ sec 




75° 

0.5 

4.5 

34.2 

27.57 

13810 

53510 

0.2440 

100° 

2.0 

8.5 

31.8 

24.94 

12500 

59750 

0.3315 

125° 

2.0 

10.0 

29.6 

24.77 

11390 

65830 

0.3846 

150° 

2.0 

9.5 

27.7 

24.63 

10410 

72060 

0.3711 


MEK <75°, 100°, 125°, 150°) 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

N 

tog 

mRo 

°C 

v*i 

sec 

cc/mln 

cm/ sec 




75° 

0.5 

4.0 

34.2 

27.57 

11060 

68860 

0.0915 

100° 

0.5 

3.0 

33.5 

27.96 

10400 

73240 

0.1520 

125° 

0.5 

3.0 

34.6 

29.57 

10240 

74380 

0.1947 

150° 

0.5 

2.5 

34.4 

29.89 

9772 

77840 

0.2439 


ETHYLENE (50°, 75°, 100°, 125°, 150°, 175°, 200°) 


Temp Semple Flow Gee Peclet 

Size Rate Velocity No* 


"to. 


mRo 


°C 


cc/mln 

cm/ sec 




VI 

o 

o 

100 

43.9 

31.88 

10110 

76980 

. 0.4247 

75° 

100 

43.6 

32.66 

9571 

81140 

0.4135 

100° 

100 

44.1 

34.18 

9239 

83910 

0.4157 

125° 

100 

45.2 

35.41 

9026 

85790 

0.4167 

o 

o 

m 

100 

44.1 

35.92 

8488 

90860 

0.3950 

175° 

100 

40.0 

33.00 

7455 

102400 

0.4683 

200° 

100 

40.0 

34.77 

7207 

105600 

0.4288 


TABLE IV 
(cont) 
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! 

ETHYLENE 100° C. Varying Sample Concentration 


Temp 

°C 

Sample 

Size 

yumoles 

Flow 

Rate 

cc/mln 

Gaa 

Velocity 

cm/eec 

Peclet 

No. 

N 

tog 

mRo 

100° 

A. 16 

143.0 

75.05 

19100 

34070 

0.3184 

100° 

4.38 

87.5 

55.45 

15400 

49060 

0.3278 

100° 

11.00 

44.7 

35.05 

9363 

82860 

0.3669 

100° 

13.30 

87.5 

55.47 

15400 

49060 

0.3297 

100° 

*** 

49.5 

37.93 

10200 

76330 

0.3160 


*** Syringe Injection of 500 








Thermodynamic Parameter, tnRo 
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Figure 11 


Effect of Temperature on Parameter 
mR / ' for Ethylene on Carbowax 1500 







PUKE COMPONENT DATA 


ON DES 


PENTANE (50°, 75°, 100°, 125°) 


Temp 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

N 

tog 

mRo 

°C 

i 

1 

sec 

cc/mln 

. _ 

cm/sec 




O 

o 

1.0 

10.0 

14.6 

2.67 

986 

638600 

0.1578 

75° 

1.0 

6.0 

21.6 

3.98 

1380 

466400 

0.2311 

100° 

0.2 

3.33 

22.9 

4.45 

1394 

462200 

0.3596 

100° 

0.5 

3.25 

22.9 

4.45 

1394 

462200 

0.3659 

100° 

1.0 

3.33 

22.9 

4.45 

1394 

462200 

0.3813 

100° 

1.0 

6.5 

27.5 

5.26 

1671 

390700 

0.3333 

125° 

0.5 

2.5 

23.2 

4.73 

1349 

476300 

0.5197 


HEPTANE (75°, 100°, 125°) 


Temp 

o_ 

Sample 

Size 

Input 

Peak 

Flow 

Rate 

Gas 

Velocity 

Peclet 

No. 

N 

tog 

mRo 

c 

1 

1 

sec 

cc/raln 

cm/ sec 




75° 

1.0 

8.0 

21. C 

3.98 

1658 

391700 

0.0431 

I -* 

o 

o 

0 

0.2 

11.0 

22.3 

4.33 

1629 

398100 

0.0806 

100° 

0.5 

11.5 

27.5 

5.26 

2005 

328600 

0.0751 

I - 1 

o 

© 

o 

0.5 

8.0 

22.6 

4.39 

1651 

393200 

0.0823 

100° 

1.0 

7.5 

27.5 

5.26 

2005 

326600 

0.0801 

100° 

1.0 

8.0 

22.6 

4.39 

1651 

393200 

0.0799 

100° 

3.0 

6.5 

27.5 

5.26 

2005 

328600 

0.0795 

125° 

1.0 

3.0 

23.0 

4.69 

1605 

403600 

0.1459 



TABLE VI 


















Figure 


2.6 


2.8 


3 


1/T x 1000 (deg K)’ 1 
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Figure 16 Effect of composition on par 
binary sample on DES at 100° 
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data on C-1500, C-102, and DES are reported in Tables VII, VIII, 
and IX. 

The activation energy of adsorption-desorption can be 
calculated from the following relation: 
mRo eC exp(-E / RT) 

which has been used to correlate data. Table X lists the acti- 
vation energies for experimental data on the three columns. 

These activation energies are in the same range as the latent 
heats of condensation for the liquid samples. This observation 
leads to the conclusion that physical adsorption is the predominant 
type of adsorption in the columns (Smith, 1970). In physical 
adsorption the activation energy for adsorption is generally 
less than 1.0 Kcal/g-mole, while the activation energy of de- 
sorption is larger because desorption is normally an endothermic 
process. This relation is shown mathematically in- this equation: 

E. « E - H 
da a 

where H , the heat of adsorption is In the same range as the 

A 

heat of condensation. Thus desorption is the controlling mechanism. 

E. Analysis of Data 

All of the experimental chromatography was done on three 
columns, the chromosorb-102 which Is packed with a mlcroporous 
polymeric bead, and the carbowax 1500 and DES columns which 
are relatively nonporous. Table XI compares the characteristics 
of the first two columns to those of the DES column, which was 
acquired this year for the purpose of separating pentane and 
heptane using a nonporous adsorbent. It is desirable to have 
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an array of columns with different characteristics so that all 
possible physical phenomena can be observed, leading to appro- 
priate general models. 

Previous workers have gathered a great deal of data for 
ethylene, pentane, heptane, and acetone on C-102 and C-1500 
(Meisch, 1973; Benoit, 1971). Comparison of these data to the 
equilibrium adsorption model simulations showed reasonable 
agreement in most cases with some obvious deviations. In general 
the data on the C-1500 column were more accurately modeled than 
those on C-102. On the C-102 column, more spreading was observed 
than was predicted by the model. It has been postulated that 
this spreading is due to pores in the packing particles not 
accounted for in the model. New data has served to clarify 
the situation. 

The largest deviation occurred when ethylene -was injected 
onto C-102 using the Carle minivalve. In this situation a 
100 yul. sample was injected as a sharp concentrated pulse with 
a duration of less than a second. The model predicted a cor- 
respondingly sharp output pulse while the actual output was a 
very wide low pulse with tremendous tailing. These results 
were shown in Figure 1. 

If porous intraparticle effects were the cause of the 
spreading then the nonporous C-1500 should show somewhat different 
behavior. When ethylene was used on this column the behavior 
was no different. A sharp input pulse was converted to a wide 
output, as is shown in Figure 18. When the concentration of 



PLOT O SIMULATED AND ACTUAL SVSTtP PE STCNSE v 0- SIMULATED OUTPUT, !• ACTUAL SYSTEM DATA 
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Figure 18 Simulated and Actual Chromatographic Data for 
Ethylene on C-1500 
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ethylene was computed it was found that the peak mole fraction 
wae nearly 0.20 In the input pulse and less than 0.04 In the 
output. The equilibrium adsorption model la not valid for large 
sample concentrations because it is simplified for small values 
of y. The exact equation 

(1 + 1/mRo) £ ylb 0 ® (1/Pe) y I b z x - 1/(1 - y) & yldz 
is simplified to a linear equation by assuming 1/(1 - y) » 1 
to yield: 

(1 + l/mRo)^y/<)G = (l/Pe) i 2 y/i) z 2 " -bylbz . 

Moreover the model is based on a linear Isotherm which may be 
a poor approximation for large concentrations. A nonlinear 
Langmuir type isotherm might better describe the situation. 

To determine the effect of sample size on this observed 
behavlor t the ethylene was Injected at lower sample pressure 
and the helium flow rate was increased. But the ethylene concen- 
tration could not be reduced enough to affect the behavior of 
the column. The Carle minlvalve can not be used to vary sample 
size over a large range. A gas tight syringe was used to get 
smaller ethylene concentrations. Figure 19 shows the results 
of this experiment. Using a wide input with low peak sample 
mole fraction) the output is extremely well modeled. The non- 
linearity appears to be composition dependent. 

Pentane and heptane on C-102 also have more spreading than 
the model predicts. When they were injected as a binary mixture 
another deviation was discovered. The value of mRo depended on 
the composition of the sample (Melsch) 1973). To discover 





whether this variable mRo was due to the porosity of the column 
packing, it was attempted to separate the mixture on the nonporous 
C-1500 for comparison. Mo separation was possible so the question 
remained unanswered. Column manufacturer's data Indicated that 
a mixture of acetone and mek would be more likely to separate 
on the carbowax column. When this mixture was tested, separation 
was minimal so It was decided to get another nonporous column 
that would separate low molecular weight hydrocarbons. 

The column chosen for these conditions was the DES column, 
a dlatomaceous earth substrate covered with di-2-ethylhexyl 
sebacate. The data from this column were extremely well modeled 
by the equilibrium adsorption model. Figure 20 compares the 
simulation to data from this column. When binary mixtures of 
pentane and heptane were used on the DES column It was found that 
linear superposition predicted the actual results quite well. 

The parameter mRo was not a strong function of composition as 
Is shown in Figure 21. Furthermore it was found that mRo was 
independent of sample size for this system. 

It was postulated that the composition dependence of mRo 
was due either to interaction between the two sample components 
or to the fact that . the. total amount of a component varied In 
samples of different composition. For example if the total 
sample was 2.0 ^1 the amount of heptane in a 0.1 heptane - 0.9 
pentane mixture would be only equivalent to about 0.2 ja.1. compared 
to 1*8 ^tl. for a 0.9 heptane • 0.1 pentane mixture; 

To clarify the behavior, pure heptane was injected onto 
C-102 with sample sizes varying from 0.2 to 3*0 ytl. 


For these 
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Figure 20 Simulated and Actual Data for Heptane 
on DES 
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runs, the value of mRo was found to be a strong function of 
sample sl 2 e (le* mole fraction in the carrier stream). Figure 
22 compares the variation of mRo with sample size for pure 
heptane runs to the variation with binary sample runs. The 
variation is similar. Any differences are probably due to the 
fact that the binary runs were made with wider input pulses , 
causing peak sample concentrations in the carrier gas to be some- 
what less than in the pure component runs of the same sample size. 
It was concluded that interactions between pentane and heptane 
had a negligible effect on mRo, and that any composition dependence 
was primarily due to the actual composition of the component in 
the gas stream. 

F. Incorporation of the Langmuir Isotherm into the Model, 

And a Solution by Orthogonal Collocation 

When the behavior of ethylene on C-102 was discussed, it 
was stated that a nonlinear isotherm mlgnt yield a more accurate 
model. This is also believed tr be the case for pentane and * 
heptane on chromosorb 102. A nonlinear Isotherm of the Langmuir 
type may possibly explain the deviations in both situations. 

The Langmuir Isotherm gives the following relation between 
the adsorbed phase and mobile phase mole fractions (Brunauer, 1943). 

X ■ Ky / (l + Ky) 

This equation is compared to the linear isotherm in Figure 23. 

For small concentrations the isotherm has nearly linear behavior. 

For large concentrations, the amount that can be adsorbed becomes 
proportionally less and finally reaches a maximum value. Larger 
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values of K make the Isotherm nonlinear at lower mobile gas phase 
concentrations. 

The nonlinearity would have the following effect on the 
behavior of the column. Large samples entering the column would 
be adsorbed proportionally less in zones of high concentration 
than in low. As a result large samples would begin to elite 
sooner than small samples and would be more spread. The quanti- 
tative effects of the isotherm on spreading and mRo can only 
be found by solving the model. 

This model takes the following form: (Appendlr Z - derivation) 

(1 ♦ l/mto(l + Ky)^)<)y/d0 «(l/Pe)<) 2 y/«) z~ - l/(l-y) 3 y/d z 

Comparing this equation to the equilibrium adsorption model, the 

2 

only difference is the term (1 + Ky) multiplying mRo. In effect, 

the parameter mRo becomes a variable dependent on composition y. 

* 

The dependence becomes stronger as K Increases. The average 

2 

value of mRo ( 1 ♦ Ky ) Increases with sample concentration. 

This is the observed behavior in experimental data. The amount 
of tailing predicted with this model can not be estimated without 
solving the equation. 

The equation must be solved numerically because it is compli- 
cated and nonlinear. Based on recent research, finite difference 
methods could be used but would probably take a prohibitive 
amount of computer time (Pfeiffer, 1972; Woodrow, 1974). The 
method of orthogonal collocation as developed by Woodrow (1974), 
has been used to advantage in some problems encountered in chrom- 
atography. This method of solution was adapted to the Langmuir 
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Isotherm problem with the some trial function and set of orthogonal 
polynomials used by Woodrow. The details concerning the formu- 
lation of the solution are In Appendix F. 

.The first set of solutions was completed for the equation 
with K set to zero* The term l/(l-y) was set equal to 1.0 and 
the equation reduced to the equilibrium adsorption model with a 
linear isotherm. Figures 24 through 26 show the results with 
/ , 7, and 15 collocation points. The solutions get progressively 
closer to the exact solution using convolution, which Is shown 
In Figure 27. However, large oscillations appear on either side 
of the main peak In the seventh and fifteenth order approximations. 

Even with 15 collocation points the main peak Is much wider and shorter 
than the exact solution. In all cases however the peak Is pre- 
dicted to occur at the proper time. 

Knowing the limitations of the problem solutions calculated 
by the method as It Is currently formulated, It was decided to 
solve the equation wltt the term l/(l-y) Included and with a 
Langmuir constant K of 100, 200, and 400. The results of these 
solutions with 15 collocation points are shown In Figures 28 
through 30. With the system parameters kopt the same as the 
other runs, these solutions are compared to the case with K equal 
to zero. As K becomes larger, corresponding to Increasing non- 
linearity in the Isotherm, the peaks elute at earlier times and 
the amount of telling appears to Increaae. Since the solutions 
e not quantitatively accurate, one can say only that they 
confirm qualitatively the expectations one has with the addition 
of the Langmuir isotherm to the model. 




PE - 9781 
MR0=0.06II 
K= 0.0 


$ 



Figure 25 Orthogonal Collocation Solution of tha Equilibrium 
Adsorption Model Using 7 Points 
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Figure 26 Orthogonal Collocation Solution of the Equilibrium 
Adsorption Model Using 15 Points 
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Figure 28 Orthogonal Collocation Solution of the Langmuir 
Isotherm Model with K « 100 











The characteristics of the collocation method as applied in 
problem were examined by performing an elgenanalysls on the 
system of ordinary differential equations resulting from the 
collocation procedure* The parameter mRo was varied from 100 
to 0*01 in the analysis and changes in eigenvalues were noted. 

The results for three collocation points are presented in Table XII. 

TABLE XII 


Elgenanalysls for Three Collocation Points 

Eigenvalue 


mRo 

Real Part 

Imaginary Part 


-1.43 

0.0 

100 

-0.883 

1.41 


-0.883 

-1.41 


-1.31 

0.0 

10 

-0.811 

1.41 


-0.811 

-1.41 


-0.724 

0.0 : 

1 

-0.446 

0.773 


-0.446 

-0.773 

0.1 

-0.132 

0.0 


-0.0811 

0.141 


-0.0811 

-0.141 

0.01 

-0.0143 

0.0 


-0.0088 

0.0153 


-0.0088 

-0.0153 


In all cases the values of the real parts of the eigenvalues 
are negative. This implies that the solutions will be stable. 

As the parameter mRo decreases however , the real part of the 
eigenvalue approaches tero which Indicates less damping of oscil- 
lations and relatively less stability. 
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If the value of mRo is set very high, the model reduces 
to the problem of pure convection, Woodrow (1974) has solved 
this equation using collocation. The solutions proved to be 
quite stable with little oscillation. Eigenanalysia predicts 
greater oscillation for low values of mRo and this was observed. 
It is reasonable to expect that a better formulation for .the 
small mRo problem exists; perhaps a different trial function or 
set of polynomials would solve the problem more accurately. 
Little work has been done in this area. 


I 
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V. CONCLUSIONS AND RECOMMENDATIONS 

From the analysis of new experimental data and observations 
documented in this report, the following conclusions are reached: 

1. Modifications to the experimental equipment have improved 
Its versatility and the accuracy of the data it provides* 

2. Changes in the data analysis programs and techniques for 
using them have added new tools to the understanding of experi- 
mental results. 

3. The thermistor probes in the chromatograph column produce 
a linear electrical response proportional to sample concentration 
at least for sample sl 2 es used in these experiments. 

A. A binary sample of pentane and heptane will not separate 
on the 100 cm chromosorb 102 column given the limitations of the 
test facility. 

5. Data from the DES column shov that samples run on this 
column are well modeled by the equilibrium adsorption model. 
Separation of the pentane - heptane binary mixture is well 
represented by the linear superposition of pure component data. 

6. Data for heptane on C-102 over a range of sample sires 
have shown that the variation of mRo with composition is a non* 
Interactive effect. A nonlinear Langmuir Isotherm is expected 
to bring the model closer to actual data. 

7. Solution of the model with a Langmuir isotherm using 
orthogonal collocation has shown qualitatively that this model 
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can predict the composition dependent mRo and the tailing found 
for samples on C-102* 

It is recommended that the Langmuir isotherm model be studied 
more closely to find out what values of K would be applicable to 
the experimental system* A method for solving the model more 
accurately should be sought* The method of orthogonal collocation 
Is extremely flexible and might be modified to better suit this 
problem* Other numerical techniques might also be studied* 



73 


VI. ACKNOWLEDGEMENT 

The author expresses his sincere gratitude to Dr. Peter 
K. Lashmet, whose patience and helpful assistance contributed 
to the completion of this project. Further acknowledgement 
Is made to Dr. G. Talne of JPL, contract monitor, and Dr. S. 
Yerazunls, project director for Rensselaer, for their Interest 
and guidance. The author acknowledges financial support from 
Rensselaer and by a graduate aselstantship under NASA grant 
NGL 33-018-091. 



74 


VII. NOMENCLATURE 

- Activation energy of adsorption 
K - Langmuir isotherm constant («■ 1/m) 
mRo - Thermodynamic Parameter 
Pe - Peclet number 
R - Gas constant 
T - Temperature 

x - Fraction of filled adsorption sites 
y - Sample concentration in the gas stream 
z « Dimensionless distance along the length of the column 
0 - Dimensionless time (vt/L) 
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IX. APPENDICES 

A. Data Processing Modifications 

1) listing of GDRP with STATMO 

2) listing of binary analysis program 

3) listing of binary simulation program MODEL 

4) listing of SIM 

B. Calculation of Mole Fractions 

C. Data Showing Thermistor Linearity 

D. Calculation of Activation ^nergy of Adsorption 

E. Derivation of Langmuir Isotherm Model 

F. Details of Orthogonal Coll >catior. Solution 

G. listing of Langmuir Isotherm Model Program 

H. Listing and Data of Eigenanalysis Program^ 


The appendices were not included in the report because of their 
length. Copies my be obtained upon request from: 


Dr. Peter K. Lashmet 
Renjselaer Polytechnic Institute 
Systems "Engineering Division 
Troy, New York 12l8l 



